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Superconductors connected to normal metallic electrodes at the nanoscale provide a potential 
source of non-locally entangled electron pairs. Such states would arise from Cooper pairs splitting 
into two electrons with opposite spins tunnelling into different leads. In an actual system the de- 
tection of these processes is hindered by the elastic transmission of individual electrons between 
the leads, yielding an opposite contribution to the non-local conductance. Here we show that elec- 
tromagnetic excitations on the superconductor can play an important role in altering the balance 
between these two processes, leading to a dominance of one upon the other depending on the spatial 
symmetry of these excitations. These findings allow to understand some intriguing recent experi- 
mental results and open the possibility to control non-local transport through a superconductor by 
an appropriate design of the experimental geometry. 



In 1964 A.F. Andreev proposed a mechanism to explain the conversion of quasiparticle currents into supercurrents at 
the interface between a normal metal (N) and a superconductor (S) [lj. In the so-called Andreev reflection mechanism 
an incident electron from the N region is reflected as a hole and a Cooper pair is created on the superconductor. This 
process can be viewed as the transfer of two electrons with opposite spins from N to S where they combine to form 
a Cooper pair. In multiterminal N/S structures a non-local version of this process where the electrons are injected 
from two spatially separated electrodes, can take place 0, 0j- The possibility to detect these crossed Andreev 
reflection (CAR) processes have attracted a lot of interest since it provides a natural mechanism to produce non- 
locally entangled electron pairs in a condensed matter device 0, @, @, HI, Q . Unfortunately, a competing mechanism 
known as elastic cotunnelling (EC) in which electrons are transmitted elastically between the two electrodes, is always 
present in actual devices yielding an opposite contribution to the non-local conductance. The theory, in fact, predicts 
that for BCS superconductors weakly coupled to non-magnetic leads the contributions of EC and CAR processes 
to the non-local conductance tend to cancel each other [lOL Spin polarisation due to magnetism in the leads 



could break this cancellation [lfj, [12|, an idea that has been explored in recent experiments [13[. Surprisingly, a 
different set of experiments by Russo et al. 14J has shown that even for the case of non-magnetic leads coupled to 
the superconductor by tunnel junctions the subgap non-local conductance can be appreciably large, exhibiting an 
intriguing behaviour in which either process can dominate depending on the applied voltage range. 

In this work we analyse the influence of electron interactions in the transport properties of this type of multiter- 
minal N/S structures. A key feature of this analysis is played by low energy collective excitations which appear in 
superconductors of reduced dimensionality 15j . When an electron tunnels into the superconductor it can excite such 
low energy modes which alters the balance between EC and CAR processes. Moreover, depending on the spatial 
symmetry of these excitations either CAR or EC processes can be suppressed, leading to a change in the sign of the 
non-local conductance as a function of the applied voltage. These findings provide an explanation to the experimental 
results of Ref. [l4| and suggest new strategies for the detection of CAR induced transport in nanostructures. 

A generic set-up for measuring the non-local resistance is shown in Fig[lji. It represents a superconducting region 
attached to three normal electrodes. Two of the leads (labelled 1 and 2 in FigOJ,) are used to inject a current while the 
voltage drop is measured on the third one. As illustrated in FigQ];, EC processes produce and injection of electrons 
into the third lead while CAR processes withdraw them from it, which yields opposite contributions to the non-local 
resistance On the other hand both processes decay on the superconducting coherence length £, which can range 
between 10 and 100 nm for typical superconductors used in experiments [13. [Til]. 

The importance of interactions in breaking the balance between EC and CAR processes can be understood by 
considering the case where the S region is sufficiently small and can be characterised by a finite charging energy E c . 
As it is sketched in FigQJ:, EC processes take place through a virtual state which will be shifted upwards by the 
Coulomb energy. The process, however, would not be blocked for any value of the applied voltage as the initial and 
final states have the same energy. In contrast, CAR processes demand that two electrons tunnel into the S region 
forming a Cooper pair, a process which requires an extra energy of AE C . Thus, the non-local conductance has a finite 
(negative) value for a voltage V smaller than 4E c /e where EC processes dominate, while it vanishes for eV > 4E C 
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FIG. 1: Top panel: a) Schematic representation of a generic multiterminal geometry where a superconducting region is coupled 
to several metallic leads and b) Double planar N/S/N junction geometry studied in Ref. [1J|- Middle panel (c): Pictorial 
description of EC and CAR processes in energy space. Lower panel (d): Feynman diagrams corresponding to the calculation of 
the non-local conductance to fourth order taking into account interactions mediated by the electromagnetic environment. Full 
lines with an arrow represent the normal and anomalous propagators and wavy lines indicate phase correlators. 



when both processes cancel each other. 

For a quantitative analysis of the influence of interactions we describe the system by a Hamiltonian H — H s + 
Hi ea ds + Ht + H env . The first three terms correspond to the electronic degrees of freedom. Hs is the usual BCS 
Hamiltonian for the S region and Hi ea d s describes the normal leads which we label with an index n. The tunnelling 
of electrons between the leads and the superconductor is described by H T = J2 n H T ^ n , with 



H T ,n=J2 MlLMin^('1^" (r1 +h.C.j, (1) 



where the integral is taken over the junction area S n , ^ (r) and ipl n a (f) are electron creation operators on the two 
sides of junction and 4> n (r) is the corresponding phase drop which is conjugate to the charge density on the junction 
(5„(r), i.e [(f> n (r), Qni^")] — ie5(r — f). The dynamics of these phase operators is determined by the Hamiltonian 
H env describing the electromagnetic environment characterising the actual experimental set-up. It is important to 
note that, due to the typical distances between the leads in the experiments, which are not much larger than £, 
correlations between voltage fluctuations on different junctions cannot be neglected, i.e. correlation functions of the 
type < 4>n4>m >, with n =/= to are non-zero. In addition, the reduced dimensions of the S region can give rise to 
the presence of collective modes within the superconducting gap, which can dominate the behaviour of the phase 
correlations. 

To obtain the transport properties of this model we use a combination of Keldysh and Nambu formalisms which is 



well adapted to analyse non-equilibrium situations in the presence of superconductivity 17[ ■ The contributions from 
EC and CAR processes to the non-local conductance G nm , i.e. the variation of the current through lead n due to a 
voltage applied on lead to, in the tunnel limit is represented by the type of diagrams shown in Fig[T]i. The solid lines 
with an arrow represent the electron propagators, while the wavy lines describe the coupling with the environment, 
i.e. they denote the phase correlators of the type < e^ n e~ % ^ m >. Let us first consider the simplest case where the 
environment can be characterised by a single electromagnetic mode of frequency luq. One can further assume that 
the leads are coupled to the S region through point contacts as illustrated in FigO Two opposite situations can be 
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FIG. 2: Pictorial representation of the effect of interactions mediated by electromagnetic modes of different symmetry on the 
non-local conductance between two point contacts. The arrows represent the phase gradient within the superconductor. While 
symmetric modes tend to suppress CAR processes the antisymmetric ones suppress the EC contribution. 



distinguished depending on the spatial symmetry of the electromagnetic mode under consideration: it can lead to 
either symmetric or antisymmetric voltage fluctuations on the two junctions. In the symmetric case, assuming that 
hu)o is much smaller than the superconducting gap A, and for S in the clean limit, at zero temperature we obtain (see 
Supplementary Information) 
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6(eV -J2 n ^o) [(-l)" 3+ " 4 cos 2 (fc F i?) - (-l)"5+»6 sin 2( fcFi? )] . 

i=l 

(2) 

Here G n ( m ) is the normal conductance of the junction n (m), R the distance separating the leads, kp is the Fermi 
wavevector, Go = 2e 2 /h is the conductance quantum and V is the voltage applied on lead m. The term proportional to 
cos 2 (kpR) corresponds to the EC contribution while the sm 2 (kpR) term arises from CAR processes. The parameter 
zq = E c /Hll>o, where E c — e 2 /2C is the charging energy on the tunnel junctions, measures the coupling to the 
electromagnetic mode. The indexes nj, with i ranging from 1 to 6 denote the number of quanta associated to the 
phase correlators on each diagram in FigQJl, 5 and 6 corresponding to the more external wavy lines. The behaviour 
predicted by Eq. @ becomes more clear in the zq <C 1 limit where it can be simplified to obtain 
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G nm ~ - " " p „ [cos 2 (fc F i?) - (1 - z 9(tuj - eV)) s\n 2 (k F R)} . (3) 
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It is worth noticing that this expression reproduces the non-interacting result [loj ] for zq = 0, where a complete 
cancellation between CAR and EC contributions takes place upon averaging over the Fermi wavelength scale. For 
finite but small zq the balance between EC and CAR is broken: for eV smaller than fuvo the CAR processes become 
suppressed and non-local transport is dominated by the EC contribution, while for eV > fnoo both contributions tend 
to cancel, as in the non-interacting case. The suppression of the CAR contribution is due to the impossibility of 
such processes to occur without producing a real excitation of the environment, as in the constant charging energy 
example. 

The situation is the opposite in the case of an antisymmetric mode. The analog of expression to Eq. @ for this case 
is obtained by an exchange of the indexes U3, with 725, ne in the last factor within brackets. As a consequence EC 
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(instead of CAR) processes will be suppressed at low voltages. The different effect of symmetric and antisymmetric 
modes is schematically illustrated in Fig[2j 

We now consider the case of a planar geometry similar to the one in the experiments of Ref. [1J] (FigjTJa) , consisting 
of a superconducting layer of thickness d > £ coupled to two normal leads by tunnel junctions. The cross section of 



the junctions in Ref. [14j was relatively large (~ 30 /ito 2 ), which allows to describe them as infinite planes |18|. The 
electromagnetic environment in this situation is characterised by the presence of propagating modes along the S/N 
junctions, which can be derived from the following model Hamiltonian 



(4) 



where the term containing the phase gradient describes the kinetic energy associated to the supercurrents in the S film, 
L being its total inductance, while the second term is the Coulomb energy of the charge accumulated on the junctions 
(assumed to be symmetric with capacitance Cp per unit area and with cross section S). In writing this Hamiltonian we 
are assuming that long range Coulomb interactions are screened by the normal electrodes acting as ground planes . 
The low energy modes which result from this model correspond to symmetric and antisymmetric voltage fluctuations 
on the junctions, with dispersion relations u>i(q) = c s \J (q tanh qd) / d and u>2{q) = c s y/q/(d tanh qd) , where q is the 
wavevector in the direction parallel to the film and c s = 1/ s/LCu- Notice that for small q the symmetric mode 
exhibits a linear dispersion with phase velocity c s while the antisymmetric one tends to a finite frequency ujq = c s /d 
in the limit q — > 0. This description of the low energy modes captures the essential features of a detailed calculation 
based on Maxwell equations for the double planar junction geometry (see Supplementary information). 

One can roughly estimate the order of magnitude of the parameters in H env for the experimental situation. Thus, 
Cp can be obtained from the typical charging energy for an oxide barrier tunnel junction E C S ~ \\ieV x fim 2 and L 
can be estimated as fJ-oX 2 /d, A being the field penetration depth [20j. The actual value of A for a Nb film is strongly 
dependent on its thickness, degree of disorder and it is also influenced by the properties of the non-superconducting 



substrate on which it is deposited [2l|. Reported values range between 100 nm and 1 \im for d ~ 10 — 100 nm |2lll22|. 
Within this range of parameters the lowest energy of the antisymmetric mode Huq can be of the order of a few meV, 
i.e. comparable to the superconducting gap in Nb, even for the smaller film thickness analysed in Ref. [14j . 

To obtain the non-local conductance Glr measured at the left interface when a voltage V is applied on the right 
junction we extend the theory developed for the single mode case, linearising with respect to the coupling parameters 
z i.2(<f) = E c /huji t 2(q), which is justified for the range of parameters estimated above. We thus obtain 



Glr = 47r ^# Ei (-2rf/0 £ (-l) a *a(S) WMq) + 1) (F( Wa (fl) + 1) + N(cj a (q)) (F(-cj a (q)) + 1)] , (5) 

F q.a=l,2 

where Gl,r is the normal conductance of each junction, F(e) = J duo^^-[f(ui + e — V) + j(uo + e + V)] is a 
thermal smearing kernel arsing from the Fermi distribution f(to) while JV(u„) is the Bose distribution function. The 
planar geometry leads to the factor Ei(-2d/£), where Ei denotes the exponential integral function. This leads to an 
exponential decay of the non-local conductance when increasing d. 

In order to understand the behaviour of Glr as a function of voltage it is convenient to first analyse the contribution 
arising from a given wavevector q. This is illustrated in the upper panel of Fig[3]for T = 10~ 2 A and Hujq — 0.4A. 
The behaviour for the different wavectors is qualitatively similar: for eV < hwi 1 hw2 the EC processes dominate, 
while CAR become more important in the voltage window Hlo\ < eV < hui2 and finally both contributions cancel for 
eV > hw2- 

The sum of all contributions yields a non-local conductance which is dominated by EC processes at V ~ * and 
decreases almost linearly with V until eV ~ Hloq- At this point the onset of antisymmetric modes produces a strong 
suppression of EC processes and there is a change of sign in Glr (CAR dominates). We can therefore associate Hujq 
with the crossing energy from EC to CAR dominated regimes. The lower panel of FigOshows the voltage dependence 
of Glr for different temperatures. As can be observed, the imbalance between EC and CAR processes driven 
by the electromagnetic modes is less pronounced for increasing temperature. The characteristic temperature for the 
suppression of the non-local conductance is set by fttoo/kB- This behaviour is in qualitative agreement with the results 
of Ref. .14?]. Moreover, the magnitude of the non-local conductance predicted by our model is in reasonable agreement 
with the experimental values. For instance, the ratio Glr/G s l r between the non-local and the direct conductances 
in the superconducting state at zero voltage and zero temperature is 0.75 Ei(— 2d/ £)(E c Su>o/hc 2 s ) which yields values 
~ 10 -3 close to the experimental ones for the parameters estimated above. A more quantitative description of the 
experimental results will be presented elsewhere. 
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FIG. 3: Upper panel: Contribution to the non-local conductance from modes with a given wavevector in the double planar 
junction geometry for feT = 0.01A and huo = 0.4A. The values of q are given in units of A/hc s - The arrows indicate the 
energy of symmetric and antisymmetric modes for q — 0.1. Lower panel: Temperature and voltage dependence of the total 
non-local conductance. The temperature values are given in units of A/fcs- The arrow indicates the energy hu>o for the lowest 
antisymmetric mode. 



In conclusion we have shown that electron interactions mediated by electromagnetic excitations lead to an imbalance 
between EC and CAR processes. Electromagnetic modes can either suppress CAR or EC processes depending on their 
spatial symmetry. Taking into account that these low energy excitations are strongly dependent on the geometrical 
characteristics of the multiterminal device, these findings open the possibility to control non-local transport processes 
through a superconductor by an appropriate design of the experimental set-up. For instance, one possibility would be 
to introduce an additional tunnel junction inside the superconducting film in the experimental geometry of Ref. [l4j ]. 
This N/S/S/N nanostructure would allow to control the dispersion relation of the electromagnetic modes by varying 
the Josephson coupling between the S layers by means of a magnetic field. Let us finally point out that the high 
sensitivity of non-local transport to the electromagnetic modes could be used as a tool to analyse these excitations in 
hybrid nanostructures. 

The authors would like to thank fruitful discussions and correspondence with D. Beckmann, A. Morpurgo, S. Russo, 
C. Urbina, D. Esteve, W. Herrera, R.C. Monreal and J.C. Cuevas. 
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